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1. Introduction 

Solid-state single photon emitters are an important step towards scalable quantum information 
technology fflHEI. a single emitter resonantly coupled to a high-finesse microcavity offers a 
natural functional unit to realize important applications in quantum cryptography and commu¬ 
nication sciences snail a. In recent years, there has been strong interest in the negatively- 
charged nitrogen-vacancy (NV ") color center in diamond for coherent control of its optical 
and spin properties |8l [91 flOl . However, although the NV center has long spin coherence times, 
of importance for many quantum applications mmmm, its fluorescence spectrum is 
stretched over a huge incoherent phonon side-band, that constitues 96% of the entire emis¬ 
sion El . More recently, other color centers in diamond were explored as a possible alternative 
to NV EUTTl; for example, the ZPL of silicon-vacancy (SiV) center (ZPL ~ 738 nm) is found 
to be much stronger than that of NV, constituting 70% of the spectrum El . Recent studies have 
also established the fundamental attributes of the SiV emitter including its electronic structure 
and polarization states 1191 [20] EU H8]| as well as single photon indistinguishability from these 
sources ll22l . validating SiV as a strong candidate for chip-based quantum optical applications. 

On the fabrication side, many important steps have been taken to integrate diamond color 
centers into chip-scale devices, including hybrid approaches |23ll24l[25ll26 :| in which diamond 
samples with color centers are juxtaposed to cavities fabricated in non-diamond materials, as 
well as monolithic all-diamond approaches [[27] [28] [29], : 30] [3T] [32] E21 in which the cavity it¬ 
self is hewn out of single crystal diamond. The latter, though technically more challenging, has 
proven to be a superior platform for diamond-based quantum photonics. On the other hand, the 
utility of third-order nonlinearity %* 3 'in diamond microphotonic devices has been explored in 
a recent work f34l which demonstrates the generation of frequency combs by parametric oscil¬ 
lators fabricated in single crystal diamond thin films. Such capabilities open up an intriguing 
possibility - to build a monolithic all-diamond emitter/frequency-converter in a single cavity 
design that efficiently collects and down-converts the color-center photon into low-loss telecom 
frequency channels for long distance communication. 

Chipscale frequency down-conversion of quantum signals has been discussed by various 
authors for different material systems f35l |36] [37l . In particular. Ref. lf37l has proposed a sili¬ 
con nitride (SiN) micro-ring resonator that can convert few-photons coherent light states from 
visible to telecom frequencies. Here, we propose a monolithic diamond structure that can be 
used to convert color-center Fock-state single photons from emission wavelengths to telecom¬ 
munication bands (~ 1.5 /im). We present a detailed theoretical description of the conversion 
process, analyzing important practical concerns such as nonlinear phase shifts and frequency 
mismatch, which were not considered in previous papers. Additionally, we present an efficient 
design technique for realizing perfect phase-matching in nonlinear optical cavities. Our analysis 
predicts sustainable power requirements for a chipscale nonlinear device with high conversion 
efficiencies achievable at total pump powers below one Watt. 

Although our theoretical analysis is generally valid for any color center, we primarily tailor 
our cavity designs to the SiV emitter. Specifically, we consider a single SiV emitter implanted 




Fig. 1. (a) Schematic of an emitter-cavity system in which a single emitter is em¬ 
bedded in a nonlinear cavity supporting four resonant modes at frequencies 
{(Boo COsc, (Otic, C°b2c}- The eigenstructure of the emitter is represented by a simplified 
three-level system with states |r), \e ), |g). A laser trigger with frequency co re and inten¬ 
sity |fl| 2 addresses the emitter states |r) and \e). The photon released from the emitter 
is collected by the cavity and down-converted to telecom through the four-wave mixing 
Bragg scattering process (FWM-BS). The latter process obeys the frequency-matching re¬ 
lation m s + a>b 2 = ft>o + ft>/,i , where <Mq and m s are the frequencies of the emitter and telecom 
signals, and ft)/,] and ft )/,2 are the frequencies of the NIR and telecom pump lasers respec¬ 
tively. (b) Alternatively, the single photon (OJo) can be down-converted by the difference 
frequency generation (DFG) process in which it is broken up into one signal (ft)j) and two 
pump photons (ft)/,], ft)/, 2 ), satisfying the frequency relation: coq — ft)/,] — ft )/,2 = ft) s . Flere, 
we choose the two pump photons to be degenerate, ft)/,/ = ffl /,2 = ft)/,. 


in a high-Q microcavity. In the case of triggered single photon emission, the electronic structure 
of the emitter can be well-approximated by a simple three-level diagram |f35l with the electronic 
states denoted by |r), \e) and |g) (Fig. QJ). To start with, the emitter is prepared in the state 
|r). An external trigger excites the system to the state \e ) which subsequently decays to the 
ground |g), emitting a photon of frequency coq into the surrounding cavity. The photon release is 
followed by frequency down-conversion (within the same cavity) via a nonlinear wave-mixing 
process which utilizes the inherent third order susceptibility ^' ^of diamond. 

A common scheme for ^ l Vl -bascd frequency down-conversion is four-wave mixing Bragg- 
scattering (FWM-BS) (36] [371, in which the incoming photon at emitter frequency coq is scat¬ 
tered by coherent pump lasers at frequencies oj/,| and ft )/,2 to yield an outgoing photon at fre¬ 
quency o s , satisfying the frequency matching condition ft) 0 + ft)/,/ = ft), + ft)/, 2 (Fig. fTJi). Al¬ 
ternatively, the down-conversion can be realized by a difference frequency generation (DFG) 
scheme in which the SiV emitter photon is directly broken up into one signal and two pump 
photons, with frequencies satisfying the relation ft)q — ft )/,2 — ft)/,] = ft), (Fig.fTJ^). In this paper, we 
present carefully designed diamond microcavities that can enhance frequency down-conversion 
via FWM-BS or DFG. In particular, we design our cavities to host four resonant modes 
(ft)0c, co sc , (Obic> COb 2 c) that respect the frequency matching conditions, (Oo c + (Ob\ c ~ ft) se + ft)/> 2 c 
for FWM-BS and ft) o c — ft)/, 1 „ — ft)/, 2 , ~ (0 SC for DFG. 

In what follows, we mainly focus our discussion on FWM-BS; however, the analysis is 
equally applicable to DFG, for which we briefly discuss design specifications, efficiency and 
power predictions. In our analysis, we take into account the effects of frequency mismatch as 
well as those of self and cross-phase modulation (ignored in the earlier work |[37lB and find that 
even in their presence, maximum efficiency can be achieved by appropriate power inputs. We 






found that the maximum efficiencies can vary from 30 % to 90%, depending on the internal 
quantum efficiency of the color center (i.e., the ratio of radiative to non-radiative decay) and 
the quality factor of the cavity mode at the emitter frequency. Our cavity designs ensure that 
these efficiencies can be achieved by input pump powers lower than one Watt. 


2. Theoretical considerations 

The effective Hamiltonian of the emitter-cavity system including the dissipative terms has the 
following form (in the Schrodinger picture): 

= ticooc djflo + ft© !C d' s a s + hcoo d ee + hco, a rr 

b (Up r ctjj h \a | 

ii,b 

+ Y J ft<*>ic{l-ttib\ab\ 2 ) a]&i 
i 

+ hCl(t)e~ ia>ret G er + hQ.* (t)e i(0ret Ore 

+ hgZVLGegdo + hgZPL&gedl ( 1 ) 

-ft/3 4aoe , ' ( ® i ’ 2 ~“ i ’ l) '-ft/T dld s e~ i ^ C0b2 ~ (0bl '> t 

+ ri Yj (glGegdi +giO ge dj S ) 



Here, d and a' are annihilation and creation operators for the quantized cavity field at the 
emitter (0) and telecom (s) wavelengths, as well as at certain other wavelengths (/) that happen 
to fall within the phonon sidebands of the emission spectrum. The modal amplitudes of the 
classical laser pumps are denoted by a /, ( b £ {b \.b2 }) and are normalized so that \ai,\ 2 is the 
energy in the b mode. (f, 7 = \i) (j\ is the atomic operator connecting the / 1 h and ;th emitter 
states. Q.(t) is the Rabi frequency coupling the trigger pump and roe transition whereas 
gZFL addresses the coupling between the ZPL decay and the 0th cavity mode. Similarly, the 
coefficients g / describe the phonon sideband transitions that coincide with the cavity resonances 
at wavelengths /. Our model also includes self and cross-phase modulation via coefficients a 
introduced by the presence of the classical pump fields, whereas the strength of frequency 
conversion is described by /3. These nonlinear coupling coefficients can be calculated from 
perturbation theory (see Appendix^- The leakage of each cavity mode is denoted by K, which 
is related to the quality factor of the mode by K = y NC denotes the spontaneous decay rate 
of the emitter into continuum loss channels, including both radiative and non-radiative decays 
|35.|. Note that in the above Hamiltonian, we have omitted nonlinear phase-modulation effects 
~ “ a 2 ll38l owing to the quantum signals themselves, as these effects are vanishingly small 

at the single-photon level. 

In addition to the Hamiltonian describing the quantum-mechanical degrees of freedom, we 
must also consider the coupled-mode equations for the classical degrees of freedom, namely 
the pump amplitudes a/, | and a^. These equations are given by [39.] : 


dab l 
dt 

dttbl 
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where P b \ and P \,2 denote the powers of the incoming lasers coupling to the cavity modes ft)/, 1 c 
and C 0 /, 2 c, with corresponding coupling lifetimes z sb i = and Z sb 2 = ■ Here, z b i = 

and T /,2 = are the overall lifetimes of the two cavity modes, which include all other sources 
of cavity losses (e.g. radiation and/or material absorption). Since one has the freedom to choose 
the frequencies (ft)/, 1 , ft)/, 2 ) of the incoming pump lasers, it is always possible to operate under 
the resonant condition for the classical pumps, 


c Obi = ®£>lc (l — — Ub\b2\ab2\ 2 ) ( 2 ) 

(0b2 = G>b2c (1 — CCbl,2\ab2\~ ~ 1 2 ) , ( 3 ) 

which allows us to ignore effects associated with self- and cross-phase modulation lf40l l39l . 
It follows that the steady-state energy in the classical cavity mode are given by \a b \\ 2 = 

Yj^-Pbi, |a& 2 1 2 = Y^Pb 2 - Linear stability analysis (40l |39l also confirms that these solutions 
are stable. 

The nonlinear coupling coefficients can be calculated from the perturbation theory ISIED 
(see also Appendix lAli: 


= T \ ^c^sc(—^-—^-)Pb\Pb2 
4 V T sb 1 T s b2 


f dv £q {Xxxyy (E* • Eq) (E/,/ • E£ 2 ) + Xxyxy (E* ■ E b \) (E() • E£ 2 ) + Xxyyx (E* ■ E£ 2 ) (Eq • E/,j)} 


a V b ~4 


VJdVe o|E ( )p y/fdVe s |E . ( | 2 y/fdV £/,i|E m | 2 V fdV £/,i|E /, 2 | 2 


3 f dV EoiXxxyy \E fl \ 2 \E b \ 2 +Xxyxy |E m -E *| 2 + Xxyyx |E^-E/,|“} 


(f dv dv e b W) 




3 fdV £ 0 {(Xxyxy + Xxyyx) |E* ■ E*|" + Xxxyy |E/,-E/,|-} 


8 


(JdVe b \E b \ 2 ) 


2\ 2 


3 fdV £ 0 {Xxxyy |E M • E * 2 | 2 + Xxyxy |E W • E w | 2 + Xxyyx |E M | 2 |E /, 2 1 2 } 
WW 4 (JdVe b i\E bl \ 2 )(JdV e b2 \E b2 \ 2 ) 


(4) 

(5) 

( 6 ) 
(7) 


The initial state of the emitter-cavity system is simply a tensor product of the emitter state 
|r) and the vacuum photonic states |Oo), |0 4 ) and |0/): 


|^(0)) = |r,Oo,0„0/). (8) 

We now proceed to solve the time evolution of under the Hamiltonian Eq. [I] It is easy to 
see that the initial state |r,Oo,O s ,0/} can only couple to states with one quantum of excitation, 
|e,Oo,O s ,0/) or \g,no,n s ,ni) , no + n s + n/ = 1. Therefore, it is sufficient to expand the total 
wavefunction as a superposition of these single-excitation states S35): 

l'E(O) = c r {t)e~ i<0rt |r,0 0 ,0 4 ,0/) 

+ c e {t)e~ lWo1 |e,Oo,0. s ,0/) 

+ co(f)«- fav \g,h,OsA) 

+ c s (t)e~ ,a>st |g,0 0 ,l s ,0/) 

+ 'Eci(t)e- i °>‘ t \g,0 0 ,O s ,\ l ). 


( 9 ) 
















Fig. 2. Schematic of angle-etched diamond ring resonator (refractive index n ss 2.41). The 
resonator has a radius R and a height h. The triangular cross-section has an etch-angle 6. 
E r components of the fundamental TEO-like modes at frequencies ( co sc , and E z com¬ 

ponents of the fundamental TMO-like modes at frequencies ( a>o c , (Obic) are also depicted 
in the picture. 


Substituting Eq. 0in the time-dependent Schrodinger equation, ih -= M' I 1 ?), leads to the 
following coupled equations of motion for the coefficients: 


c r = —i£l*(t)c e 

Ce = -i£l(t)c r - igZPLCO - ij^gi c i 

i 



Co = -idoc-CQ — igZPL.Ce ~ >P*C S -— Co 

C s = -i5 sc c s -il5c 0 -^c s 

Kl 
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5oc = (Dqc (l — c(obi\ a bi\ 2 — otowkwl 2 ) — too 
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( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 


The nonlinear detunings Sq c and 8 SC , as defined by Eqs. fl5T4T6l include the nonlinear phase 
shifts introduced by the pumps and are to be distinguished from the bare cavity detunings 
coot- — (Oo and co sc + (Obi c ~ (Ooc ~ (Obic without the amplitude-dependent phase shifts. 

In the regime where the trigger intensity ~ Q. and trigger rate ~ ^ are sufficiently small 
compared to the decay rates of the system, the probability coefficients c,-, i € {e,0,,v,/}, follow 
the dynamics of c r . Thus we can adiabatically eliminate the e.Q.s.l degrees of freedom (see 
Appendix [0, and arrive at the following expression for the conversion efficiency (defined as 
the percentage of triggers for which the excitations are ultimately down-converted to telecom 

ESI): 


f dt K ss | | 2 

Jdt (RA || 2 + kb|c 0 | 2 + L/ k'/lc/l 2 + y NC |c e | 2 ) 

_ 4C\P\ 2 KqK ss _ 

161/31 4 + 4\P | 2 ((2 + C)k q k s - 8 SoeSsc) + (4 8 2 c + (1 + C) k 2 ) (45 2 + k 2 ) ' 


(17) 

(18) 


Here, C is the well-known cooperativity factor: 

















which describes the relative strength of photon capture into the cavity Oth mode over the com¬ 
peting loss channels. Eq.[l8]can be maximized by setting = 0, = 0, which yields: 

<20) 

under the constraints, 

<5oc rit = (4VT+C\W-(l+C)K QK j. (21) 

A special situation arises when 


<S erit 

°0c- 


= <5 C 


= 0 , 


( 22 ) 


in which case fi is constrained by. 


^(i+c^vlgr (23) 

This scenario was already discussed in Ref. If35l . More generally, we can algebraically 
solve Eq. |2T] together with Eqs. [2|-[3] to yield a set of critical powers and frequencies 
{ft)™ 1 , P™ 1 , ft)™ 1 , Pjp lt }- However, it must be noted that the above algebraic equations some¬ 
times lead to negative or imaginary answers. Such cases indicate that the maximal efficiency is 
no longer given by Eq.[20l instead the critical efficiencies and powers are numerically found by 
directly setting 


dF dF 

- - = 0 , -- 

dP b i dP b2 


(24) 


In the following sections, we will perform these calculations in the context of carefully designed 
realistic cavity systems. 


3. Diamond Micro-ring Resonator 

In this section, we consider concrete and realistic cavity designs for single photon frequency 
conversion based on a triangular ring resonator. We choose the triangular cross-section to take 
advantage of our recently developed angle-etched technique on monolithic bulk diamond l33l . 
The main advantage of angle-etched technique is its consistent yield and scalability compared 
to diamond-on-insulator thin-film techniques. 

The most challenging task in our cavity design is to satisfy the frequency-matching relation 
CUoc + (Obic ~ + ft)/, 2 C . One cannot straightforwardly apply conventional dispersion engineer¬ 

ing techniques developed for FWM with closely-spaced frequencies l42l since, in our case, the 
conversion spans two very different frequency bands, NIR and telecom. One can do an ex¬ 
tensive numerical search in the design space lf37ll but such an approach is time-consuming 
and computationally very expensive. In contrast, in the Appendix 0 we describe a numeri¬ 
cal/visualization technique that greatly simplifies the process of designing ring resonators that 
satisfy the frequency-matching condition. One such cavity is shown in Fig. [2] The ring has a ra¬ 
dius R = 14.4 fim, a height h = 414 nm and a slant angle 9 = 59°. Out of its many resonances, 
we will focus our attention on two fundamental transverse-electric-like (TEo-like) modes in the 
telecom band, (m s = 95,A SC = 1.613 /im) and ( m b i = 79, K\c = 1 .806 fim), and two funda¬ 
mental transverse-magnetic-like (TMo-like) modes in the NIR band, (mo = 262. Ao,- = 738 nm) 
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Fig. 3. Density plot of conversion efficiency F (defined in Eg. 1181 over the pump powers 
Pb\ and P b 2 for the cavity system discussed in Section [3] Efficiency contours are overlaid 
on the plot for easy visualization. They help identify the regime of pump powers necessary 
for high-efficiency conversion. 

and (nib 2 - 246. A/, 2 , = 776 nm), see Fig. |2] Finite-difference time-domain (FDTD) simula¬ 
tions reveal that the four modes have ultra-high radiative quality factors in excess of 10 8 . 
However, in actual experiments, the overall Q’s will be limited by fabrication imperfections 
or coupling losses. Therefore, we will examine the maximum efficiencies for two different rep¬ 
resentative cases, Q ~ 10 s and Q ~ 10 6 . For the sake of simplicity, we will assume that the 
overall Q’s are coupling-limited, so that Q, « (f v ,. i £ {0.,v./;l .b2). If this is not the case, the 
maximal efficiency will suffer by a factor of Q s /Qss (see Eq. 120b . The third-order nonlinear 
susceptibility of diamond (Kerr coefficient) is X'^ = 2.5 x 10' 21 ^4 143 M . The nonlinear phase- 
shifts a and coupling coefficient are found to be: 



Oom =0.114 (fiTT 1 , ao i2 =0.425 (/tJ) _1 , 
Osbi =0.215 (juJ) -1 , a sb2 = 0.126 
a bb i = 0.099 (juJ)' 1 , a hh2 = 0.208 (juJ)' 1 , 
«MM = 0.114 (jUJ )- 1 . 


(26) 

(27) 

(28) 
(29) 


The maximum acheivable efficiency (Eq. l20l i is determined by the cooperativity factor C 
(Eq.UD, which can be recast as: 



( 30 ) 












detuning A 0c = (u;o c - w 0 ) /uj 0c 



detuning A sc — 4 (w sc + cq, 2 c — w 0 c — w&ic) / (w sc + u\, 2 C + woe + w i>ic) 

Fig. 4. (a) Critical pump power (solid lines, left-axis) and critical efficiency (dashed lines, 
right-axis) vs. bare-cavity fractional detuning Ao r = 1 — a\)/(OQ c . Pump power plotted here 
is the total pump power, f’FjLj = +P™ 1 . Two regimes can be identified: a regime with 

multiple critical points (light red and dark blue colors) and a regime with only one critial 
point (dark blue only). The former corresponds to the regime where Eq. [20] is valid and 
the detuning can be compensated by critical pump powers. The latter corresponds to the 
regime where Eq. [20] is no longer valid and the efficiency falls off rapidly with detuning. 
(See also text.) (b) same as above except for the detuning A JC = y 4 ■ 


where 7 is the total spontaneous decay rate of the color center in the isotropic diamond medium 












4p2 

in the absence of the cavity. The quantity Y,i ~yy re P resen ts the relative strength of radiative 
decay into phonon sidebands in the presence of the ring resonator as opposed to the total decay 
rate in the absence of the resonator. Generally, one can compute this quantity by simulating the 
decay rates of all of the ring resonator modes falling within the phonon sidebands. However, in 
the case of SiV with a relatively weak phonon sideband, we can expect much of the sideband 

4g~- 

to be suppressed in between the ring resonator modes so that L / Fy 1 ■ Therefore, it follows 
that 

or (31) 

(32) 

■yNC 4„2 

since y- < 1 in general. In our calculations, we will simply take the lower bound C = . 

Assuming that the SiV is located at the electric-field maximum and that the dipole element is 
aligned with the field polarization, the coupling strength is given by gzpL = ,Uy|,| | Enlax . Sponta- 

neous emission rate J/\ y \, is related to the dipole element via ' m °de volume 

of the Oth cavity mode is Vq c = / = 252 (. The branching ratio of the SiV ZPL 

£( r maxJ|^max| \ n J 

line is » 70% US). A variable factor in our calculations is the internal quantum efficiency 

of the SiV emission (or the ratio y™ 1 //); the current literature suggests it could vary wildly from 
sample to sample, from as low as 9% l44l to as high as near 100% at low temperatures ETI . 
probably depending on material quality and method of fabrication. In Table. [7] we have com¬ 
puted maximal efficiencies and critical pump powers (Eq. [20] and Eqs. 12214231) for different 
internal quantum efficiencies and cavity lifetimes. A maximal conversion efficiency of 87% in 
the best possible scenario where y 1 ' ad /y = 1 and Q = 10 6 is limited by a moderate cooperativity 
C « 211 (see Eq. Ea.[20l>. However, we note that this limit only reflects the current experimental 
constraints on diamond fabrication (Q ~ 10 6 ); simulations have indicated that radiative qual¬ 
ity factors can be in excess of 10 7 , in which case the maximum efficiency quickly approaches 
100 %. 
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j? max 

10 s 

0.09 

2 

0.110 

0.424 

0.26 

10 s 

1 

21 

0.255 

0.507 

0.65 

10 6 

0.09 

19 

3 x 10~ 4 

0.036 

0.63 

10 b 

1 

211 

0.001 

0.037 

0.87 


Table 1. Critical powers and efficiencies. 

Next, we examine the scenario with Q = 10 6 and y- = 1 as the best possible performance 
for the given design. Figure [3] shows a density plot of the conversion efficiency (Eq. ITSl) as 
a function of the two pump powers. As observed in the figure, there exist narrow strips of 
high-efficiency regime (F > 0.85) around the pump power co-ordinates (1 mW, 35 mW) and 
(15 mW, 12 mW). A closer examination reveals two additional critical points (in addition to 
the one given in Table. [7} with efficiencies E c " 1,1 = 0.87, F cnt ’ 2 = 0.82 and pump powers 
P“ iU « 14 mW, P™ 1 ’ 2 « 11 mW and P™ u « 5 mW, P™ 1 ’ 2 « 28 mW respectively. In con¬ 
trast to the tabulated critical point, these extra points have a small but non-vanishing phase 
mismatch 8 ^ 0 and are obtained from a complete solution of the algebraic equation (Eq. 1211 


















or from the direct optimization over pump powers (Eq. [24}. The existence of multiple critical 
points is a result of a complex interplay between two competing decay channels, (i) the decay 
of the excited state with the concomitant photon release into mode 0 and (ii) the conversion of 
the released photon to the telecom mode s. If the leakage from mode 0 to s is too strong, no 
cavity enhanced emission occurs, degrading the efficiency, whereas if this leakage is too weak, 
nonlinear conversion efficiency vanishes. These decay rates are controlled by the nonlinear de¬ 
tunings, <5() c , 5 vc, and nonlinear coupling j3 parameters, which are in turn controlled by the input 
powers. Specifically, Sq c tends to inhibit the emitter —> 0 decay channel, 8 SC tends to inhibit the 
0 —> s decay channel, and /I tends to enhance the 0 —> s decay channel. The complex interplay 
between these parameters introduces multiple critical points over the I\ i —Pya plane. 

As we have shown above, the nonlinear detunings <5’s are directly controlled by incident 
pump powers, which can, therefore, be appropriately chosen to ensure maximum conversion 
efficiency. The magnitudes of the critical pump powers in turn depend on how much frequency 
detuning the system begins with. Theoretically, the cavity parameters can be numerically op¬ 
timized so that these detunings are made as small as possible. However, random errors in 
the actual fabrication process can introduce significant deviations from any designated cav¬ 
ity frequencies. Therefore, it is advisable to consider the effects of frequency detunings due to 
possible variations in ring radius R and height h. For this purpose, we define the /?«re-cavity 
fractional detunings: A 0c - = 1 - ftto/®0o A* c = (o^+o^+o^+t^l)/4 ■ Figure 0] shows the total 
critical power n™' tal = PS 1 ' + and corresponding critical efficiencies F cnt as a function of 
A( )( . (Fig.[4ja)) and A JC (Fig.[4jh)). (Note that F cnt denotes all critical points, including local and 
global maxima.) Here, A( lr is varied by varying (Oq while keeping all other frequencies fixed, 
and A„. is varied by varying (O sc . We find that generally there exist two regimes of operation: 
a regime where there are two critical efficiencies, one at F cnt = 0.87 (light red dashed line) 
and the other at F cnt < 0.87 (dark blue dashed line), and a regime where only the smaller effi¬ 
ciency survives, F mt < 0.87 (dark blue dashed line). These two regimes clearly reflect to what 
extent the detunings can be compensated for by incident pump powers. In the first regime, the 
detunings can be completely compensated and the maximal efficiency is only limited by coop- 
erativity C. The global maximum is explicitly computed from Eq. [2T] as F tnl = 0.87 and can 
be realized by applying either of the two critical powers (light red solid lines) corresponding to 
either vanishing or non-vanishing 5’s. The required power then grows with the increasing A’s 
while the global maximal efficiency stays the same. In the second regime, Eq.[2T]breaks down 
(leading to complex solutions), and the critical parameters can only be computed from setting 
Tffr = 0’ 77777 = 0- Clearly, in this regime, the detunings cannot be fully compensated so that 
the maximal efficiency falls off rapidly. 

4. Difference frequency generation 

So far, we have tailored our resonator design to the FWM-BS process. Designing a cavity for 
the alternative scheme of DFG is significantly more challenging since the pump modes must be 
designed at a wavelength of ~ 3urn, more than thrice that of the SiV emitter, so that the gen¬ 
erated signal falls within the telecom band (see Fig.0. For these wavelengths, phase-matching 
can only be achieved by utilizing higher-order waveguide modes; as a consequence, the overlap 
between the cavity modes suffers, degrading the /3 coefficient. Nevertheless, with the help of 
the phase-matching method described in Annendixlcl we were able to obtain realistic designs 
that still yield high efficiencies at reasonable pump powers. Such a design is shown in Fig. 0 

In order to separate TE-like and TM-like bands as wide as possible, we choose a shallow etch 
angle of 9 = 70°; the radius and height are 34 u m and 463 nm respectively. The waveguide di¬ 
mensions are chosen to be larger than the previous design for FWM-BS in order to host high-Q 
modes at very different wavelengths: (mo = 461, = 0.738 jUm), (m s = 235, = 1.567 jx m) 
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Fig. 5. Left panel: a diagrammatic representation of difference frequency generation pro¬ 
cess (DFG), Right panel: E r and E z components of a higher-order TE-like mode (for the 
emitter photon), a fundamental TMO-like mode (for the signal photon) and a fundamental 
TEO-like mode (for the pump photons). TE and TM characters of the modes are ill-defined 
since the mirror symmetry in z (the out-of-plane direction perpendicular to the plane of the 
resonator) is strongly broken. Each mode possesses appreciable in-plane and out-of-plane 
electric field components, which lead to a non-vanishing /3 (see Eq.[4]>. 

and (mi,\ = mf ,2 = mb = 113, i c = A^c = A& = 2.789 /dm). (Note the pumps are degenerate 
so that we are actually dealing with three cavity modes instead of four.) Here, it should be noted 
that even though we make use of two TE-like modes and one TM-like modes (see Fig.0, we 
found that the overlap does not vanish because TE and TM characters of the modes are ill- 
defined for the triangular waveguide, all the modes having siginficant in-plane and out-of-plane 
components. From FDTD, radiation Q’s are found to be in excess of 10 8 ; we assume opera¬ 
tional Q’s of 10 6 . The mode volume of the Oth cavity mode is Vo, « 1052 ^ j , yielding a 
cooperativity factor of about 50 (assuming 100% internal quantum efficiency of the emitter). 
The analysis of DFG exactly follows that of DFWM, except that we interchange CObi — (O/a 
and Eb\ While the self- and cross-phase modulation coefficients are comparable to those 

in the case of FWM-BS, the overlap coefficient is about 2.5 x 10~ 4 (/rJ) -1 , which is more than 
two orders of magnitude smaller than in the case of FWM-BS. With these numbers in place, we 
calculated the critical efficiency and powers which are found to be 64% and ~ 1 W. 

5. Conclusion 

We have presented efficiency and power requirements for frequency down-conversion of SiV 
single photons to telecom wavelengths. Our cavity designs are tailored for a potential all¬ 
diamond monolithic approach which directly utilizes the nonlinear optical properties of di¬ 
amond. This stands in contrast to more common hybrid approaches 051 which employs an 
external nonlinear material, usually a X' 21 medium, for frequency conversion. Although the 
power requirements (on the order of 10 — 100 mW) are greater for all-diamond systems since 
the inherent /( ''nonlinearity of diamond is weaker, the maximum-achievable efficiencies are 
comparable to ^ l2 *-based systems, for example Ref. 051 . The advantage of a potential all¬ 
diamond approach is the superior quality of the native emitters in monolithic diamond as well 
as relative ease of fabrication, scalability and high throughput promised by the advent of angle- 
etched techniques j45l . 







A. Nonlinear terms in Hamiltonian 

The perturbative nonlinear interaction energy is given by j46l 

SU NL = -^J dV E*(co)-P NL (ft)). (33) 

There are detailed theoretical methods on how to rigorously quantize a nonlinear system 1 381 
l46l but we will not follow those arguments here since we are only interested in the perturbative 
regime. 

The nonlinear polarization for four-wave mixing Bragg scattering process is given in ll43l . 

P™ L (co s ) = Deoxjj ] kl (co s ; co 0 , ,-(o b2 )E 0j E blk E; 2l , (34) 

where D is the permutation factor and in this case D = 6 11431 . The electric fields shoulqd be 
normalized as 


f dV = 

(35) 

J dV Ej((o) • eE fe (co) = )j\a b \ 2 

(36) 


In an isotropic medium like diamond, the components of the nonlinear susceptibility tensor 
are given by Xijki = XxxyySijSki + Xxyxy$ik$ji + Xxyyx$u$jk- Substituting Eq.[34]in Eq.[33]leads 
to the expression for [3 as given in Eq. Q] Similar calculations can be performed to yield the 
expressions for a’s, Eqs. [5]{7] For wavelengths of interest which are well into non-resonant 
regime, full-permutation symmetry holds ll43l and we approximate that Xxxyy ~ Xxyxy ~ Xxyyx ~ 
Xxxxx/3 = Z (3) /3. For diamond, X {3) = 2.5 x 1(T 21 m 2 /V 2 |43). 

B. Adiabatic elimination 

When fl is sufficiently small, c r becomes the slowest dynamical variable in the system leading 
the other three. In this case, c e , cq . c s , c/ all decay approximately at the same rate as c r . In 
other words, c e , co, c s , ci are simply proportional to c r (they follow c r except for a constant 
pre-factor). To obtain these proportionalities, we can formally set c e . co, c s , c; = 0 individually 
(adiabatic elimination) and solve for these variables as functions of c r . This leads to a simple 
rate equation for c r : 


Cr = (i8r- y)c r , (37) 

where 8, is some phase factor giving rise to unitary oscillations and K, is the effective rate of 
population loss predicted from the state |r). The latter is given by 

Kr - l^l 2 (\ _ 4\P\ 2 Ck 0 k s _ 

^ \ 16|/3| 2 +8|^| 2 (koK^(1+ C) — 45o c 5 ic ) + ^45g c + (1 +C) 2 (45^. + iCj) 

Self-consistency of adiabatic elimination requires that the population loss out of r does not 
exceed the rate K s with which a photon can leak out through the cavity mode s, K, < k s . Ad¬ 
ditionally, in our calculations, we have rigorously checked the validity of adiabatic elimination 
by direct comparison with exact numerical integration of Eqs. IT() 1 - [T61 using |Q| ~ 0.1 tec,. 




C. Ring resonator design 

We begin by examining the modes of an infinite straight waveguide having a triangular cross- 
section. Specifically, we chose an etch-angle of 60° and a height of 0.25a where a is an arbitrary 
normalization length to be chosen later. The modes of the waveguides can be quickly computed 
by standard band structure solvers; for our simulations, we use the freely-available open-source 
MPB PP71 which compute eigenfrequencies / (in units of c/a) at specified k\ (in units of 2 ft/a). 
We select two bands which are far-separated in frequency but still have significant overlap; such 
bands are readily afforded by the fundamental TEO and fundamental TMO modes (see the inset 
of Fig. | 6 j. Next, we define the phase mismatch as a function of frequencies (for the case of 
FWM-BS) 


Skifi ,/ 2l /3,/4) = h +k 2 -k 3 - k 4 . (39) 

We want two of the frequencies, say, f 2 and f 4 to lie within the NIR band and the other two, f\ 
and / 3 , within the telecom band. Therefore, we identify: k\ = kTEo(/t)i k 2 = &tmo(/2 ): k 3 = 
^TEO {f3 )} k 4 = £tmo(/4)- 

8 k is a function of four independent variables f \, f 2l / 3 , f 4 which is hard to visualize. We 
subject the frequencies to the constraint fx+f 2 = / 3 +f 4 ■ In order to reduce the number of free 
variables further, we make / 3 a function of f\ and f 2 , i.e., / 3 = / 3 (fx,f 2 ). This also makes f 4 
a function of f\ and f 2 . The functional form of / 3 can be chosen arbitrarily. A simple intuitive 
choice would be to set = fx — 8 f for some offset Sf. Since the number of variables has 
been reduced to two, we are ready to visualize 8 k. Choosing df = 0.1, we plotted a contour of 
8 k = 0 over the f\ vs. f 2 plane (see Fig. [ 6 ]). Each point along the contour yields four waveguide 
modes with frequencies that are perfectly phase-matched for an efficient FWM-BS process. In 
particular, we chose a point on the contour (red circle in Fig. [ 6 j: fx =0.95, f 2 = 2.00, / 3 = 
0.85, f 4 = 2.10. It is easily noticed that f 4 and f\ have the approximate ratio as two frequencies 
from SiV and telecom bands. So we identify f 4 ‘- = /o = TrfWJim' yi^Tug a = 1.55 jinx and 
f s =f\ | = l 63 /i m ■ This also specifies the height of the waveguide, h ss 390 nm. Remembering 
that ring resonances satisfy the condition m, = InRn^ jf,. i G { 0 .,v./?l . 82 ), we can find an 
optimal R which makes m, an integer (or closest to an integer) for each of the four modes. 
We found that R ~ 14.7/rm approximately gives mo = 262, rn s = 95, m/,i = 79, 1111,2 = 246, 
also satisfying the condition mo + m/,\ = m s + mi >2 . Now we can proceed to simulate the ring 
resonator via full-scale FDTD, using the predicted design parameters. A quick inexpensive 
search (on the order of ten iterations) in the vicinity of the above parameters directly yields the 
design given in Section^ 

A similar procedure is followed to design a cavity for DFG except that the angle is made 
shallower 0 = 70°. We found that the 8 th order TE-like band together with fundamental TE- 
like and TM-like bands yield a regime of phase-matched frequencies (not shown) similar to the 
one depicted in Fig. [ 6 ] thereby identifying appropriate parameters which lead to the resonator 
design of Section|4] after a quick numerical search. 





fa/c 

Fig. 6. Phase-matching diagram for the FWM-BS process in a triangular waveguide of 
h = 0.25 a and 6 = 60°. The red circle indicates the frequencies used in the design of ring 
resonator. Inset: fundamental TEO and TMO bands computed by MPB. a is an arbitrary 
normalization length to be chosen later (see text). 









